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We present a model where sterile neutrinos with rest masses in the range ~ keV to ~ MeV 
can be the dark matter and be consistent with all laboratory, cosmological, large-scale structure, 
as well as x-ray constraints. These sterile neutrinos are assumed to freeze out of thermal and 
chemical equilibrium with matter and radiation in the very early Universe, prior to an epoch of 
prodigious entropy generation (“dilution”) from out-of-equilibrium decay of heavy particles. In this 
work, we consider heavy, entropy-producing particles in the TeV to ~ EeV rest-mass range, 
possibly associated with new physics at high-energy scales. The process of dilution can give the 
sterile neutrinos the appropriate relic densities, but it also alters their energy spectra so that they 
could act like cold dark matter, despite relatively low rest masses as compared to conventional dark 
matter candidates. Moreover, since the model does not rely on active-sterile mixing for producing 
the relic density, the mixing angles can be small enough to evade current x-ray or lifetime constraints. 
Nevertheless, we discuss how future x-ray observations, future lepton number constraints, and future 
observations and sophisticated simulations of large-scale structure could, in conjunction, provide 
evidence for this model and/or constrain and probe its parameters. 

PACS numbers: 95.35.-|-d; 14.60.Pq; 14.60.St; 95.30.Cq 


I. INTRODUCTION 

In this paper we show how sterile neutrinos with rest 
masses in the ~ keV to ^ MeV range could evade all 
existing cosmological and laboratory bounds, comprise 
all or a significant component of the dark matter, and 
behave as cold dark matter (CDM). The idea of an elec- 
troweak singlet (sterile neutrino) as dark matter is not 
new [1-38]. Some of these models run afoul of x-ray ob¬ 
servations or large-scale structure considerations or both, 
as discussed in Refs. [3, 39-42] and also in Sec. IV below. 
However, many of them still remain viable. 

Most of these models posit no sterile neutrinos at ex¬ 
tremely high temperatures in the early Universe and en¬ 
gineer a build-up of sterile neutrino density via active 
neutrino scattering-induced decoherence, lepton number- 
driven medium enhancement of that process, or particle 
decay. In a different class of models [10, 11, 20, 22, 31], a 
population of thermally decoupled sterile neutrinos forms 
in the early Universe with a density below the equilib¬ 
rium density, and the subsequent expansion of the Uni¬ 
verse reduces both the density and the momenta of sterile 
neutrinos, making them acceptable dark matter candi¬ 
dates. It is also possible that scattering-induced deco¬ 
herence produces a population of sterile neutrinos with 
a number density that is initially too high, but is subse¬ 
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quently reduced to an acceptable level by entropy gen¬ 
eration (“dilution”) effected by out-of-equilibrium decay 
of heavy sterile neutrinos different from the dark matter 
candidate sterile neutrino [9]. 

Finally, one can start with the dark matter candidate 
sterile neutrinos in thermal and chemical equilibrium at 
some high temperature [23, 24, 26, 30]. These sterile neu¬ 
trinos eventually decouple, and, after decoupling, have 
their number density reduced to the required dark mat¬ 
ter density by dilution due to out-of-equlibrium heavy 
particle decay. Our paper follows this approach. 

In this paper we consider heavy, unstable dilution gen¬ 
erators, henceforth referred to as “dilutons,” with rest 
masses in the ^ TeV to ^ EeV range. There is noth¬ 
ing sacred about this diluton mass range. It could be 
higher and it could be somewhat lower, though cosmo¬ 
logical considerations discussed below may limit how low. 
In any case, our model requires large dilution and there¬ 
fore, the diluton lifetime against decay has to be long 
enough for the dilutons to survive until the Universe has 
cooled to temperatures well below their rest mass. 

A diluton particle with those properties could be, for 
example, another sterile neutrino, with an extremely 
small vacuum mixing with active neutrino flavors. The 
small vacuum mixing would enable evasion of laboratory 
neutrino mass and accelerator bounds, and cause a rela¬ 
tively long lifetime [43, 44]. This sterile neutrino diluton 
perhaps could be one of the heavy right-handed neutrino 
species invoked in the seesaw mechanism for explaining 
neutrino mass phenomenology. 

Another possibility is that the diluton is a supersym¬ 
metric particle that decays into standard-model particles 
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via R-parity-violating couplings. The relatively long life¬ 
time required for this particle could be effected by hav¬ 
ing R-parity be a nearly respected, but ultimately broken 
symmetry [45]. However, this would also mean that dark 
matter would not be a lightest supersymmetric particle 
(LSP). A supersymmetric diluton which is an LSP would 
require that the scale of supersymmetry be very high. 

As mentioned above, the idea of diluted-equilibrium 
sterile neutrino dark matter (DESNDM) has been dis¬ 
cussed before for particular kinds of dilutions. For exam¬ 
ple, gauge extensions of the standard model containing 
right-handed neutrinos, i.e., left-right symmetric mod¬ 
els, can furnish diluton candidates decaying around the 
weak decoupling epoch (temperature T ~ 1 MeV) [23], 
or the QCD scale (T ~ 100 MeV) [24]. These mod¬ 
els can produce ~ keV rest-mass-scale sterile neutrino 
warm dark matter. Variants of this model posit more 
massive dilutons, decaying above the electroweak scale 
(T ~ 100 GeV), and can produce heavier, colder dark 
matter [26, 30]. 

In this work we shall remain agnostic about the iden¬ 
tity and rest mass of the diluton, with the only assump¬ 
tions being that it also has an equilibrium distribution 
prior to its decoupling and subsequent decay, and that 
all of its final decay products thermalize. This allows us 
to make sterile neutrinos in the rest-mass range ~ keV 
to ^ MeV be CDM. 

Entropy generation in the early Universe can have con- 
strainable consequences. In particular, it has been shown 
that entropy injection at or near the weak decoupling 
scale (T ~ MeV) can be constrained by big bang nu¬ 
cleosynthesis (BBN) and radiation energy density (Aeff) 
considerations [44, 46-49]. However, there would be no 
effect on BBN and A^eff if particle decay-generated en¬ 
tropy injection occurs sufficiently prior to weak decou¬ 
pling, so long as all of the diluton decay products ther¬ 
malize in the plasma of the early Universe. Evading other 
potential cosmological bounds may argue for an even 
higher temperature scale for a significant dilution event. 
For example, though we do not know where the baryon 
number is made, it could be produced at or above the 
electroweak scale. Our dilution event may have to occur 
above the temperature epoch associated with baryogene- 
sis; otherwise, a higher pre-dilution baryon number would 
have to be produced, placing potentially unattainable de¬ 
mands on baryon generation mechanisms. Nevertheless, 
we also consider cases where the dilution event occurs 
at temperatures much lower than the electroweak scale, 
and these would require an appropriate accompanying 
baryogenesis scheme. 

In Sec. H, we briefly discuss various mechanisms for 
sterile neutrino production and thermalization and de¬ 
scribe how dilution can be incorporated into the thermal 
history of the Universe. This is followed by an assessment 
of the effects of dilution in Sec. HI. Section IV gives an 
overview of the various possible means to observation- 
ally or experimentally probe our model and constrain its 
parameters. We conclude in Sec. V. 


II. STERILE NEUTRINOS AND THE HISTORY 
OF THE EARLY UNIVERSE 


Figure 1 illustrates the thermal history of the early 
Universe and highlights key epochs in the history of ac¬ 
tive and sterile neutrino species. This figure gives a par¬ 
ticular example of when dilution events might occur in 
viable scenarios of diluted equilibrium sterile neutrino 
dark matter (DESNDM). The dilution event is depicted 
in this figure to occur before the electroweak phase tran¬ 
sition, which might be more favorable from the point of 
view of some baryogenesis models. However, we also con¬ 
sider longer-lived dilutons in this paper. 
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FIG. 1. Cartoon illustrating how the DESNDM model fits 
together with the thermal history of the Universe. So long 
as the diluton Dh decays exclusively into standard model 
particles, and the process goes to completion before weak de¬ 
coupling, the decay products can completely thermalize. In 
such a scenario, all standard model particle spectra remain 
thermal, and neither big bang nucleosynthesis nor Aefr are 
affected in any way. However, since the sterile neutrinos Vs 
are already decoupled, their number densities are diluted rela¬ 
tive to the particles still in equilibrium. The active neutrinos 
eventually decouple, and subsequently, the pairs annihi¬ 
late away as the Universe cools, further diluting the active 
and sterile neutrino number densities relative to the photons. 
Scenarios with dilution built in prior to the electroweak scale, 
such as the one shown here, are likely better suited to meeting 
baryogenesis requirements. 
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A. Sterile neutrino production in the early 
Universe 


Sterile neutrinos can interact subweakly with ordinary 
matter via their mixing with active neutrinos [50]. For 
example, in a simplified model where the sterile neutrino 
mixes with only one of the active species, the interaction 
rate of a sterile neutrino with the background plasma in 
the early Universe can be estimated as 

( 1 ) 


where Gp and T are the Fermi coupling constant and 
the plasma temperature, respectively, and 9m is the effec¬ 
tive active-sterile in-medium mixing angle. At high tem¬ 
peratures, these in-medium mixing angles can be heav¬ 
ily suppressed due to active neutrino interactions with 
background matter. At a plasma temperature T, the ef¬ 
fective in-medium mixing angle for a neutrino state with 
momentum p satisfies 


sin^ 29m 


(p) sin^ 29V 

A2(p) sin^ 29y + [A(p) cos20„ — Ud — Ur]^ ’ 


( 2 ) 

where 9y is the active-sterile mixing angle in vacuum, 
and A(p) = Am^/2p, with Am^ being the appropri¬ 
ate mass-squared splitting in vacuum. Vd and Vr are 
the finite-density and finite-temperature potentials felt 
by the active neutrino. In the absence of an appreciable 
net lepton number, the finite-density potential is negligi¬ 
ble. The finite-temperature potential, for T < Mw (Mw 
is the lU-boson mass), is given by Ur = —G^gpT"^, where 
Geff can be taken to represent some overall neutrino in¬ 
teraction strength, summed over all the particle species 
in the background plasma. 

In the limit of negligible lepton number, i.e., Ud ^ 0, 
the thermal potential Ur causes heavy suppression of the 
effective in-medium mixing angle at high temperatures. 
Numerically, the thermal potential Ur and the vacuum 
oscillation term A(p) can be calculated to be 


(3) 


For a keV-MeV rest-mass-scale sterile neutrino, jUrj ~ 
A(p) at T ~ O.I-l GeV. At higher temperatures, the 
thermal potential dominates, and the scattering rate goes 
like Fj/^ oc whereas at lower temperatures, where 

the vacuum oscillation term dominates, F^^ oc T^. The 
sterile neutrino scattering rate is therefore maximal in 
these intermediate temperature regimes, where the vac¬ 
uum and thermal terms are comparable in magnitude. 
For equilibration, this rate has to be greater than the ex¬ 
pansion rate, which assuming radiation-dominated con¬ 
ditions is given hy H = 5 ^/^ T^/rUpi. Here 

g is the statistical weight in relativistic particles and the 
Planck mass is nipi. With these considerations it can be 


shown that sterile neutrinos with vacuum mixing angles 
smaller than sin^ 20„ ~ 10“® (10 keV/mg) can never be 
in thermal equilibrium with the plasma, so long as their 
interactions arise solely through these mixings. There¬ 
fore, in such scenarios, their relic density cannot be set 
via a freeze-out process, unlike ordinary active neutrinos. 

Such mixing with active neutrinos can, however, lead 
to sterile neutrinos being produced athermally in the 
early Universe via scattering-induced decoherence, as was 
first pointed out by Dodelson and Widrow [1]. Ac¬ 
tive neutrino scattering-induced decoherence into ster¬ 
ile states is driven by the active neutrino scattering rate 
F^^ ^ G|. T^. However, this scattering also gives rise to 
the active neutrino matter potentials described above, as 
well as quantum damping, both of which serve to inhibit 
active-sterile neutrino conversion. The sterile neutrino 
production rate via this mechanism can also be shown to 
peak at temperatures of ^ 0.1-1 GeV. 

Combined constraints from x-ray and Lyman-a ob¬ 
servations rule out a Dodelson-Widrow-type sterile neu¬ 
trino being all of the dark matter (although it may com¬ 
prise some fraction of the total Hdm, depending on the 
mass and mixing angle). If the lepton number is sizable, 
however, the scattering-induced conversion rate can be 
resonantly enhanced [2, 3]. Sterile neutrinos which are 
resonantly produced in this way can have considerably 
lower vacuum mixing angles for a given relic density and 
can have “colder” energy spectra, as compared to the 
Dodelson-Widrow case, and can therefore evade some of 
these bounds. The Lyman-alpha bounds are also relaxed 
if dilution takes place after the production of sterile neu¬ 
trinos through this mechanism [9]. 

Alternatively, one could envision ways to bring the 
sterile neutrino into equilibrium early on by invoking 
some additional interactions, such as, e.g., left-right sym¬ 
metric models, suitably broken at some energy scale, 
above which a right-handed Z-boson effects equilibration 
of this sterile neutrino with the plasma [23, 24, 26, 30]. 
We adopt this approach in our paper, assuming that the 
dark matter candidate sterile neutrino attains an equilib¬ 
rium distribution via some beyond-standard-model inter¬ 
actions, before freezing out. Additionally, for simplicity, 
we also assume the diluton to be thermally populated, 
either through similar or through different interactions 
as the dark matter candidate sterile neutrino. 

There are a number of models that would produce an 
equilibrium population of sterile neutrinos through pro¬ 
cesses other than oscillations. In the context of grand 
unified theories (GUT) with the popular S'O(IO) gauge 
group, the right-handed neutrino transforms as a compo¬ 
nent of a 16-dimensional representation that also includes 
all the standard model fermions. At the GUT scale, the 
right-handed neutrinos can be produced in equilibrium 
by the exchange of the GUT scale bosons. These inter¬ 
actions freeze out at lower temperatures, and the sterile 
neutrinos go out of equilibrium, but their distribution 
functions remain the same and scale with temperature. 
Depending on the mode of S'O(IO) breaking, there may 
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or may not be any entropy production between the GUT 
scale and the electroweak scale. 

As an example, if the breaking occurs along the route 
S'O(IO) —>■ SU{5) X U{1) with a subsequent breaking of 
SU{5) down to the standard model at a scale close to the 
S'O(IO) breaking scale, the population of sterile neutrinos 
can exist with near-thermal distribution functions. In the 
case of an alternative symmetry-breaking route, leading 
to a left-right symmetric SU{2)l x S'[/(2)/{X ... model, 
the right-handed SU (2)^^ gauge bosons can keep the ster¬ 
ile neutrinos in equilibrium down to the scale a few orders 
of magnitude below the SU{2)r breaking scale, which 
has to be at the TeV-scale or higher [51]. If the right- 
handed neutrinos are doublets of SU{2)r, the Majorana 
mass cannot be greater than the SU{2)r breaking scale; 
hence, we “naturally” get a Majorana mass much smaller 
than the Planck scale. Furthermore, one might expect to 
find the right-handed Z and W at the Large Hadron Col¬ 
lider (LHC) [52]. 

More generally, any high-energy theory that includes a 
gauge U{1)b-l either as a subgroup of the GUT group 
(e.g., S'O(IO)), or as a stand-alone feature (e.g., split see¬ 
saw model [22]) generates an approximately thermal pop¬ 
ulation of sterile neutrinos through the exchange of the 
U{1)b-l gauge boson. This is because our sterile neutri¬ 
nos carry a lepton number, and hence a nonzero B—L. In 
the case of a split seesaw model of Ref. [22], a first-order 
phase transition was employed to dilute the density of 
sterile neutrinos, but such dilution can be small or none 
if the U{1)b-l breaking transition is not strongly first 
order. 

In all of these scenarios with additional interactions 
at high temperatures, equilibration of a massive ster¬ 
ile neutrino species will lead to their “overproduction” 
in the early Universe, if the sterile mass is greater than 
rus ~ 100 eV. In that case, entropy generation after the 
sterile neutrinos have decoupled can help dilute their relic 
densities to a level consistent with astrophysical data. 


B. Thermal decoupling 


The DESNDM model assumes that at very early times, 
the keV-MeV mass sterile neutrino, as well as the dilu- 
ton are in thermal and chemical equilibrium with matter 
and radiation. A particle species thermally decouples 
from this plasma when its scattering rate falls below the 
expansion rate of the Universe (or equivalently, its mean 
free path becomes longer than the Hubble length). This 
could happen, for example, if the interactions responsi¬ 
ble for keeping these particles in equilibrium in the very 
early Universe weaken considerably, or cease to operate 
below a certain energy scale, e.g., following the breaking 
of some symmetry as in the above examples. 


C. Decay-induced dilution 


In our model, we assume that the heavy diluton de¬ 
cays into standard model particles (such as photons, pi- 
ons and a variety of charged and neutral leptons) after 
both the diluton and the dark matter sterile neutrino 
have decoupled, thus pumping vast amounts of entropy 
into the plasma. Let us define S = s ■ as the co¬ 
moving entropy of the plasma, where s = (27r^/45) gs 
is the entropy density. Here, (/„, a and T are the effec¬ 
tive entropic degrees of freedom, the scale factor, and the 
plasma temperature, respectively. The symmetries of a 
Friedmann-Lemaitre-Robertson-Walker (FLRW) space- 
time imply that the comoving entropy is conserved as 
long as there are not any timelike heat flows. However, 
out-of-equilibrium particle decays can source such heat 
flows. The rate of change of comoving entropy resulting 
from diluton decay is given by 


dt 


nn ■ a mn , 

rji JT ^ 

TH ± 


( 5 ) 


where niH, th, and nn are the mass, lifetime, and 
number density, respectively, of the diluton, and /t is 
the fraction of the total mass energy of decay products 
which thermalizes in the plasma. Here, we are assuming 
T <C niH, so that the decaying particle’s energy can be 
approximated as its rest mass. Assuming that the dilu- 
tons were initially in thermal equilibrium, their number 
density relative to photons is given by 


where we have assumed a relativistic Fermi-Dirac shaped 
energy/momentum distribution function for the diluton 
with a temperature parameter Th, and with g = 2 spin 
degrees of freedom. This assumption about the shape of 
the diluton distribution function is tantamount to an as¬ 
sumption that the diluton particles are relativistic when 
they decouple. Putting together Eqs. (5) and (6), we can 
write 

\ds 135 c(3) ^ Oils' -t/rn .71 

S dt 47r4 T ^ 

where we have used S/{a^n-y) = tt^ gs/(45C(3)). This 
allows us to numerically compute the entropy added to 
the plasma at each time-step, for a given diluton mass 
and lifetime. Assuming that the diluton decays exclu¬ 
sively into standard-model particles, and that the decay 
occurs well before active neutrino (i.e., weak) decoupling, 
all the decay products can fully thermalize (i.e., /r = I), 
thus preventing some of the entropy from “leaking away” 
into decoupled particles. Defining the “dilution factor” 
F = S'finai/'S'initiai, as the ratio of the comoving entropies 
of the plasma after and before the dilution event, we can 
write 


gsa^T^ = gs,ia^T^F. 


( 8 ) 
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Here, the subscript “i” (initial) appearing on the right- 
hand side of the equation is supposed to indicate the 
onset of dilution, whereas the quantities on the left-hand 
side (without this subscript) are to represent any post¬ 
dilution epoch. The dillution factor can be > 10, so 
that most of the entropy of the Universe is generated in 
this decay-induced dilution scenario. The DM candidate 
steriles, which are assumed to decouple prior to diluton 
decay, do not benefit from any of this entropy injection, 
and the temperature parameter that characterizes 
their energy distribution simply redshifts inversely with 
the scale factor, i.e., a = Ui Assuming = Ti 
(i.e., steriles initially in equilibrium), we can write 



Temperature [MeV] 


T \gs,^ ' Fj 


(9) 


The “cooling” of the sterile neutrino sea relative to 
the plasma can, therefore, be seen to be a combined ef¬ 
fect of dilution from particle decay, as well as the dis¬ 
appearance of statistical degrees of freedom as the Uni¬ 
verse cools. The latter effect is most prominent across 
the quark-hadron transition at T Ri 170 MeV, where gs 
drops sharply by a factor of ^ 3. Overall, gs decreases 
from ^ 100 at T > 100 GeV, to gs ~ 10.75 by the time 
the active neutrinos decouple at T ~ MeV, at which point 
the only relativistic degress of freedom are the photons, 
electrons, positrons, neutrinos and anti-neutrinos. Sub¬ 
sequently, as the temperature drops significantly below 
the electron rest mass, the pairs annihilate away trans¬ 
ferring most of their entropy to photons, which offsets 
the temperature of active and sterile neutrinos relative 
to the photons by a factor of r; (4/11)^/^. As a result, 
gs ~ 43/11 at late times, taking into account the entropy 
in photons, neutrinos, and anti-neutrinos. 

Since the number densities of sterile neutrinos and 
thermally coupled particles are ~ Tj^ and ~ T^, respec¬ 
tively, this process results in an effective decrease in the 
number of sterile neutrinos relative to the plasma, di¬ 
luting their relic density down from an initial thermal 
value. And as long as the diluton decay happens suffi¬ 
ciently early and all decay products completely thermal- 
ize in the plasma, no imprint is left on Ngg or BBN. 

Figure 2 illustrates the process of dilution by plot¬ 
ting the active and sterile neutrino temperature curves 
against the plasma temperature. Figure 2a (top) gives 
examples of dilution events occuring prior to the elec- 
troweak scale, whereas Fig. 2b (bottom) considers post- 
electroweak dillution cases. The sterile neutrino mass nis 
is calculated in each case so as to get the appropriate relic 
density for them to be the dark matter, as explained in 
the following section. In particular, the case of mg Rs 7.1 
keV has attracted recent interest in light of certain x-ray 
observations (see Sec. IV B). However, the range of appli¬ 
cability of this model extends beyond just that particular 
case. 



Temperature [MeV] 


FIG. 2. Active neutrino (solid, red) and sterile neutrino 
(dashed, green and dot-dashed, blue) cooling curves as a func¬ 
tion of plasma temperature, illustrating how a combination of 
out-of-equilibrium particle decay and loss of statistical degrees 
of freedom can dilute decoupled particles relative to species 
still in thermal equilibrium. The curves begin to separate as 
the diluton starts decaying, pumping entropy into the plasma 
and diluting the the number density of the sterile neutrino 
sea. The figures are for the cases with diluton rest mass rriH 
and lifetime th, and associated dilution factor F, and dark 
matter sterile neutrino rest mass m^, as labeled. The sterile 
neutrino mass, for a given diluton mass and lifetime, is cho¬ 
sen to give closure fraction Ds = 0.258 for Hubble parameter 
h = 0.6781 (in units of lOOkms”^ Mpc“^). 


III. CONSEQUENCES OF DILUTION 


A. Dark matter particle mass and relic density 


We have seen how dilution from out-of-equilibrium par¬ 
ticle decay subsequent to sterile neutrino decoupling, but 
prior to active neutrino decoupling, will dilute the dark 
matter steriles relative to the active neutrinos. In turn, 
e^-pair annihilation subsequent to active neutrino decou¬ 
pling cools the active and sterile neutrinos relative to the 
photons by the usual factor of Ri (4/11)^/^, so that the 
ratio of sterile neutrino temperature to photon tempera- 
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ture Tj at much later epochs is 


k 

where 5s,wd ~ 10.75 is the statistical weight for entropy 
at the active neutrino decoupling epoch. 

Since the sterile neutrinos would be expected to be 
nonrelativistic at the present epoch, their energy density 
would simply be the product of their number density and 
rest mass, ms- However, since these particles 

would decouple while they were still relativisitc, their en¬ 
ergy distribution would retain its relativistic Fermi-Dirac 
shape, with a temperature parameter that redshifts 
inversely with the scale factor. We can, therefore, write 
p,, = [(3C(3)T^J/(27r2)] - ms, where C(3) « 1.20206 
is the zeta function of argument 3, and where we add 
in both right- and left-handed sterile states. This im¬ 
plies that the sterile neutrino rest-mass contribution to 
closure is fig = Pi^Jpcrit, where the critical density is 
Pcrit = 3 Hq mpj/STT, and where Hq and rupi are the Hub¬ 
ble parameter at the current epoch, and the Planck mass, 
respectively. 

Using Eq. (10), and given the observed cosmic mi¬ 
crowave background (CMB) temperature, T^o = 2.725 K, 
and a Hubble parameter h = Hq/ (lOOkm s“^ Mpc“^), 
the sterile neutrino rest mass which would account for a 
closure parameter at the current epoch is 


d 5s,wd 1 
11 gs,t F 


n J./0 


( 10 ) 


■nis 


11^ Wpiffg 5,,, 
16C(3) T^q 5s, wd 


2.26 keV 


9s,^|gs ,wd 

10 




( 11 ) 




FIG. 3. Figures showing contours of the sterile neutrino rest 
mass ms (in keV) which would give a current relic abundance 
of Us = 0.258 for Hubble parameter Hq = 67.81 km/s Mpc“^ 
in the DESNDM model, plotted against a parameter space 
spanned by diluton rest mass mn and decay lifetime th- 


Figure 3 shows contours of sterile neutrino rest mass 
(in keV) that would give the measured dark matter relic 
abundance, for different diluton rest masses and lifetimes. 
Figure 3a (top) explores diluton rest masses in the PeV- 
EeV range, whereas Fig. 3b (bottom) has dilutons with 
TeV-PeV rest masses, but longer lifetimes. 

In summary, the diluton rest mass and lifetime to¬ 
gether determine the dilution factor (i.e., the ratio of 
final to initial comoving entropy), which in turn picks 
out an appropriate sterile neutrino rest mass for them to 
be the dark matter. 


B. Dark matter collisionless damping scale 

Dark matter particles can be classified as “hot,” 
“warm,” or “cold,” depending on their collisionless damp¬ 
ing scale. For a sterile neutrino distribution characterized 
by a temperature parameter , the collisionless damp¬ 
ing (free-streaming) length scale, comoving to the current 


epoch, can be estimated as 


Afs ~ 


X 


0.27 Mpc 



7 -I- In 


( ^ 


T^\ ({PITU\ 

T ) \ 3.15 ) 

3.15 0.14 Y 

■ (p/T),, ■ nraK^)\ ’ 


( 12 ) 


where we have adopted the analysis from Ref. [53], with 
slight modifications. Here, 17^ is the fraction of en¬ 
ergy density at the current epoch contributed by all 
nonrelativistic matter (dark matter -|- baryons). In the 
DESNDM model, the average ratio of momentum to tem¬ 
perature for the decoupled sterile neutrino energy distri¬ 
bution is given by the thermal value, {p/T)^^ « 3.15, 
since the spectrum retains its thermal, Fermi-Dirac shape 
postdilution, albeit with a significantly reduced temper¬ 
ature parameter. 

Although some of the sterile neutrino rest masses dis¬ 
cussed here may appear low enough to be in trouble with 
some of the more stringent bounds based on observed 
structure in the Lyman-a forest [54], this is not the case. 
The reason is that these diluted sterile neutrinos would 
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have small collisionless damping lengths, since the tem¬ 
perature parameter that characterizes their energy distri¬ 
bution is lower than the photon temperature. Therefore, 
these sterile neutrinos effectively behave as if they were 
thermal particles with a higher mass, = nis (T 
so that 


Afs 


0.27 Mpc 



7-fin 




Here, the superscript “cd,” short for “collisionless 
damping,” is being used to emphasize that this thermally 
adjusted particle mass is the relevant parameter that de¬ 
termines the collisionless damping scale. Using Eqs. (10) 
and (11), this can be expressed as 


m 


cd 

s 


= rus 



9s,t 
9s,wd 


1/3 


^1/3 


a; 18.5 keV 


9s,i/gs ,wd 

To 


4 

3 




(14) 


For all but the lightest sterile neutrinos in the rest- 
mass range under consideration, this effective mass is 
above the ~ ISkeV limit where damping of structure 
from dark matter particle free streaming could come into 
conflict with observation [54]. Therefore, in this model, a 
relatively light (rUs a few keV) sterile neutrino can also 
function as cold dark matter. This is because the high 
dilution factor that is required to get the correct relic 
density, also serves to suppress the free-streaming scale. 
The total mass inside the sterile neutrino free streaming 
scale, Mps = (47r/3) Apg pm, can be calculated as 


Mps 


4 X 10® Mq 


10 -f In 


/20keVy h'^\ 
V mf ) V 0.14 ) 
mf 0.14 A1^ 
20keV ’ rimh^ )\ ’ 


(15) 


where Pm is the total energy density in nonrelativistic 
matter (baryons + dark matter) at the present epoch, 
corresponding to closure fraction Mps defines a 

mass scale for fluctuations, below which they would ex¬ 
perience considerable damping via dark matter particle 
free streaming. Again, it must be emphasised that the 
effective thermally-adjusted particle mass that sets 
this scale is not the sterile neutrino rest mass, but is 
scaled relative to it by the ratio of photon temperature 
to sterile neutrino temperature. 

Figure 4 depicts how the total mass inside the sterile 
neutrino free-streaming scale varies as a function of ster¬ 
ile neutrino rest mass, in this model. In Fig. 4a (top), 
we consider the case where the diluted-equilibrium ster¬ 
ile neutrinos are all of the dark matter, i.e., Qjgh? = 
Udm ~ 0.12, whereas in Fig. 4b (bottom), we al¬ 
low the closure density parameter fig h? to vary, in or¬ 
der to account for situations where these sterile neutri¬ 
nos may not be all of the dark matter. Aside from the 


relatively unimportant logarithmic dependence, the vari¬ 
ation of A/ps with iTis and Ug /i^ can be quantified as 
Mfs oc mj^ ■ Hs h^. 


Thermally adjusted effective mass (keV) 




FIG. 4. Top: total mass inside the sterile neutrino free- 
streaming (i.e., collisionless damping) scale, Mfs, as a func¬ 
tion of sterile neutrino rest mass rris, for fig — 0.12. Also 
shown are the temperature ratio /T.y (on the y2 axis) as a 
function of ms, for fls h? = 0.12, and the corresponding ther¬ 
mally adjusted effective mass (on the x2 axis). Bottom: 
contours of Mfs, in solar masses, plotted against sterile neu¬ 
trino rest mass ms on the x axis, and sterile neutrino closure 
density parameter Qsh^ on the y axis. 


For a broad range of sterile neutrino rest mass and di- 
luton properties, the DFSNDM model can produce what 
is effectively CDM, at least as far as the existence of 
a Lyman-a forest is concerned. However, in certtain 
ranges, the sterile neutrino character may not be entirely 
CDM-like, from the point of view of certain other aspects 
of structure-formation. For example, for sterile neutrino 
rest masses in the range ^ 5-10 keV in the DFSNDM 
model, collisionless damping scales can be ^ 10 ^Mq. 
These could fall in a range of interest for the core/cusp 
problem and other issues in dwarf galaxy morphology 
under current investigation. Further discussion follows 
in Sec. IV D. 






















C. Matter-dominated epochs in the early Universe 

The presence of a heavy particle with a number den¬ 
sity comparable to that of thermal particles raises the 
possibility of an epoch of matter domination in the early 
Universe. For example, a ^ PeV rest-mass-scale diluton 
with a lifetime of order 10“^^ s will linger around until 
the temperature of the plasma has dropped to about a 
100 GeV. This means that there will be a period of time 
where the total energy density of the Universe is dom¬ 
inated by the diluton rest mass. This is illustrated in 
Fig. 5 (dot-dashed, blue curve). 

10 ® 

10^ 

^ 10^ 

^ 10“ 

a 

^ 10 ^“ 

I 10-‘ 

S 10-“ 

S jq-8 
10-10 
10 -'“ 

10“ 10' 10“ 10“ 10^ 10“ 10“ 10' 10“ 

Temperature [MeV] 

FIG. 5. Curves showing the total horizon mass energy (solid, 
red), the Jeans mass (dashed, green), and the ratio of the 
energy density contributed by diluton rest mass to the total 
energy density in radiation (dot-dashed, blue, plotted on y2 
axis), as a function of plasma temperature. The Jeans mass 
can be seen to drop relative to the horizon mass in the matter- 
dominated epochs. The features (bends) in the Mj and Mhor 
curves at T 100 MeV are a consequence of the relatively 
sudden change in relativistic degrees of freedom across the 
QCD transition. The diluton rest mass and lifetime in this 
example are niH = 2.91 PeV and th ~ 0.7406 x 10“^^ s, re¬ 
spectively. 

Because of this effect, the causal horizon proper length 
starts to increase, from a value of (ihor(i) = 2t in 
radiation-dominated conditions, towards (ihor(0 = 3< in 
matter-dominated conditions. In a radiation-dominated 
Universe, the total mass energy contained in the causal 
horizon, Mhor = (dir/S) ptot, is a factor of few smaller 
than the total Jeans mass, Mj = (Jtt/S) Aj ptot, which 
is the mass scale above which gravitation can overcome 
pressure support and cause fluctuations to grow in am¬ 
plitude. Here, ptot is the total energy density, and the 
Jeans length is Aj = Cg rupi/^ptot, with Cg = l/v^ as 
the sound speed (even when the energy density is dom¬ 
inated by the diluton rest mass, it is the thermally cou¬ 
pled, relativistic particles in the plasma that provide the 
pressure support, and hence determine the sound speed). 

An epoch of early matter domination can cause the 
Jeans mass to drop below the causal horizon mass scale, 
meaning perturbations that are in causal contact can 


now start to grow in amplitude. For our test case 
of a 2.91 PeV rest-mass diluton with a lifetime th = 
0.7406 X 10“^^ s, this phase lasts for less than a decade 
in temperature (see Fig. 5), i.e., a few Hubble times. 
For this particular choice of parameter values, the rel¬ 
evant mass scale of the fluctuations in this regime is 
1O“^-1O“®M0. Given the small horizon mass scale, the 
limited interval of matter domination, and significant ra¬ 
diation content of the plasma, it is unlikely that any non¬ 
linear regime fluctuations produced in this epoch can sur¬ 
vive to later epochs with appreciable and constrainable 
amplitudes [55-58]. 


IV. OBSERVATIONAL AND EXPERIMENTAL 
HANDLES 

Broadly speaking, the parameter space of sterile neu¬ 
trino rest mass and vacuum mixing angle can be con¬ 
strained by dark matter stability considerations, vari¬ 
ous x-ray observations, and also through kinematic argu¬ 
ments, e.g., phase space considerations [59], and bounds 
on dark matter collisionless damping from Ly-a forest 
observations. Some of these constraints are summarized 
in Fig. 6 and described in further detail in the subsequent 
subsections. 


A. Dark matter stability considerations 

For any particle to be an acceptable dark matter can¬ 
didate, it has to be stable against decay or annihilation 
processes over the lifetime of the Universe (to 4 x 10^^ 
s). Massive sterile neutrinos are able to decay by virtue of 
the posited mixing with active species in vacuum. For a 
sterile neutrino in the keV-MeV range, the predominant 
disappearance channel is a tree-level decay mediated by 
a Z boson, into three active neutrinos, i.e., Vg —>■ Su. For 
a sterile neutrino that mixes with a single active species, 
the decay rate is given by 

sin^ ( 16 ) 

Requiring that the lifetime Tg « > 10^® s 

leads to the following constraint on the sterile neutrino 
mass-mixing parameter space 



If the sterile neutrino rest mass is greater than twice 
the electron rest mass, i.e., mg > MeV, then that opens 
up an additional decay channel Vg —>■ ve'^e~. The rate for 
this channel is one-third the decay rate into iv, resulting 
in a factor of 4/3 enhancement in the overall decay rate, 
making the above constraint slightly more stringent at 
these relatively higher rest masses. 
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FIG. 6 . Sterile neutrino rest mass and vacuum mixing pa¬ 
rameter space, as constrained by x-ray/ 7 -ray observations 
(constraints from Refs. [60-65] shown here in purple), Ly- 
a forest limits on collisionless damping (from Ref. [54], shown 
here in dark red), as well as the requirements that dark mat¬ 
ter is stable (black) and does not rethermalize after dilu¬ 
tion (gray). Also shown are the regions of allowed param¬ 
eter space that can produce some observable effects on dwarf- 
galaxy morphology (green, cross-hatched region corresponds 
to Mps ~ 1 O®- 1 O®M 0 in the DESNDM model), and on pulsar 
kicks (sky-blue, diagonally hatched region reproduced from 
Ref. [21]), respectively. The x-ray/ 7 -ray constraints shown 
here are premised on sterile neutrinos being all of the dark 
matter. The rest-mass ranges encompassed by the Ly-Q and 
dwarf galaxy regions shown here are specific to the DESNDM 
model. Finally, the gray solid lines are contours of sterile 
neutrino closure fraction produced by the Dodelson-Widrow 
mechanism: Ddw = 0.26 (right) and 0.0055 (left). 


B. X-ray observations 


In addition to the aforementioned decays via tree-level 
weak processes, there is a radiative electromagnetic decay 
branch arising via one-loop interactions [3, 39, 43] which 
provides a photon with energy ms/2. The decay rate for 
this electromagnetic branch is given by 


r 




sin^ 20„ a Gp ^ 
6.8 X 10"^^ s“^ 


9TOs \ 


20487r4 ) 


/ sin^ 26y \ 

/ TOs \ 

1 10-10 j 

V 1 keV / 


(18) 


Although this rate is 0{a) smaller compared to the 3i/ 
channel, the fact that it leaves an electromagnetic imprint 
makes it much more viable for indirect detection [39, 40]. 
Consequently, this radiative decay channel has been used 
to give the currently most stringent constraints on many 
models for sterile neutrino dark matter. X-ray obser¬ 
vations of the Milky Way [60-62, 66 , 67], Andromeda 
(M31) and other local group galaxies [63, 68-70], dwarf 
spheroidals [71-74], and galaxy clusters [64, 75-77], as 
well as measurements of the diffuse and unresolved cos¬ 
mic x-ray backgrounds [65, 78], have been used to con¬ 


strain the parameter space of rest mass and vacuum mix¬ 
ing for sterile-neutrino dark matter models. 

The model proposed here can evade these constraints, 
again effectively because of the dilution involved in cre¬ 
ating their relic densities. Since the active-sterile mixing 
angle has no bearing on the relic density in this model, 
it can be made arbitrarily small. Since the decay-photon 
emissivity from the sterile neutrinos is proportional to 
the square of the appropriate vacuum mixing angle, all 
the above bounds could therefore be evaded. 


1. The 3.55 keV x-ray line 


Recent analysis [79-81] of x-ray emission from various 
sources has led to the detection of a previously uniden¬ 
tified monochromatic x-ray emission at a photon energy 
of around 3.55 keV, possibly arising from electromag¬ 
netic decay of a 7.1 keV rest-mass sterile neutrino into 
an active neutrino and a photon [82]. From the observed 
fluxes, and with the assumption that the sterile neutrinos 
constitute all of the dark matter, the inferred best-fit vac¬ 
uum mixing angles are sin^ 20„ Ri 7 x 10“^^. While the 
existence of this line in terms of statistical significance, as 
well as its interpretation as having a dark matter origin 
are still up for debate [83-94], the possibility remains in¬ 
triguing, and various sterile neutrino dark matter models 
can have their parameters tailored to fit this particular 
scenario [29, 30, 32-36, 38]. Future x-ray telescopes such 
as ASTRO-H and ATHENA, as well as microcalorimeter 
sounding rocket experiments such as Micro-X [95, 96], 
with their high energy resolution, could help settle the 
verdict on this case one way or the other [97]. 

Some of the results presented in this paper, such as 
in Fig. 2, as well as the discussion in Secs. IIIC and 
IV C, have used this posited 7.1 keV rest-mass sterile 
neutrino as an example. However, much of the analysis is 
also applicable more generally, for a wide range of sterile 
neutrino parameters. 


2. Looking for heavier sterile neutrino dark matter 


We have shown that our model can dilute an initial 
thermal distribution down to the right relic density even 
for much heavier sterile neutrinos, whose electromagnetic 
decay branches would fall outside the purview of tele¬ 
scopes such as Chandra, XMM-Newton or Suzaku. How¬ 
ever, some of this higher rest-mass range would lie in 
a suitable regime for other x-ray/ 7 -ray telescopes, e.g., 
Fermi-GBM, which can probe the 20-50 keV rest-mass 
range for sterile neutrinos [62], or NuSTAR, which is de¬ 
signed to see x-rays in the 3-79 keV range [77, 98], cor¬ 
responding to TOs = 6-158 keV, or INTEGRAL, looking 
at 18 keV -8 MeV photon energies [61, 67]. 
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C. Dependence on mixing angle 

In our model, the dark matter relic density is set by 
assuming an equilibrium distribution of sterile neutrinos 
in the early Universe, followed by an epoch of out-of- 
equilibrium heavy particle decay, which engineers an ap¬ 
propriate amount of dilution. Therefore, our model does 
not rely on the active-sterile mixing angle, as far as set¬ 
ting the relic density is concerned. However, in certain 
regimes, the mixing angle can have other important con¬ 
sequences, and can therefore be constrained. 

1. Avoiding re-thermalization after dilution 

The DESNDM model rides on the assumption that the 
sterile neutrino decouples in the early Universe, followed 
by an epoch of entropy injection which cools the sterile 
sea relative to the plasma. However, if the sterile neutri¬ 
nos were to re-thermalize after this, the entire purpose of 
dilution would be lost. 

The sterile neutrino scattering rate in the plasma 
would have to be greater than the expansion rate at 
some epoch for them to re-thermalize, and as demon¬ 
strated in Sec. HA, this would not happen for small 
enough active-sterile vacuum mixing angles, as long as 
the lepton number is negligible. This puts an upper limit 
on the active-sterile vacuum mixing angle of sin^ 29y < 
10“® (10 keV/rUs), for nis ^ keV-MeV, if the dilution 
model is to work. In practice, vacuum mixing angles in 
this range are already ruled out by x-ray observations 
and/or constraints on sterile neutrino lifetime, so the up¬ 
per limits arrived at using thermalization arguments are 
merely of academic interest. 

2. Warm dark matter component produced by 
scattering-induced decoherence 

As mentioned in Sec. HA, sterile neutrinos can be pro¬ 
duced in the early Universe via scattering induced deco¬ 
herence, even in the absence of a lepton number. So long 
as the sterile neutrinos mix with the active neutrinos, this 
process is unavoidable, and in the context of our model, 
it can produce an additional contribution to the sterile 
neutrino relic density. Also, if the dilution event were to 
happen prior to T ~ 0.1-1 GeV, the sterile neutrinos pro¬ 
duced by scattering-induced decoherence would possess a 
higher average kinetic energy compared to their diluted- 
thermal-relic counterparts. This decoherently produced 
component will therefore contribute a warm tail to the 
overall sterile neutrino energy spectrum, leading to an 
increase in the effective dark matter collisionless damp¬ 
ing scale, with likely implications for structure formation 
models. 

For a sterile neutrino with mass rris and vacuum mixing 
angle 9y with the active neutrinos, the contribution of 
this component to the closure density of the Universe in 


the zero-lepton number limit is given by [1, 21] 



where ‘DW’ is an acronym for Dodelson-Widrow. 
Clearly, one must have Hdw < f^DM ~ 0.26, and this 
puts an upper limit on the mixing angle as a function 
of sterile neutrino rest mass, in order to avoid overabun¬ 
dance of these steriles. For the particular case of the 
nis = 7.1 keV sterile neutrino described in Sec. IV B, 
Fq. (19) becomes 

f^DW-0.94(^^^). (20) 

In scenarios where the 7.1 keV sterile neutrino is all 
of the dark matter, i.e., Hs = Hdmj the inferred mix¬ 
ing angle from the observed flux of the x-ray line is 
sin^ 29y « 7 X 10“^^. This implies sin^ « 1.75 x 10“^^, 
which using Eq. (20) gives Hdw ~ 0.0055 « 0.021 0dm, 
for 0 dm = 0.26. Thus, if the 7.1 keV sterile hinted at 
by x-ray observations were to be all of the dark mat¬ 
ter, then about 2% of its total number density (in the 
zero-lepton number limit) would be produced by the 
Dodelson-Widrow mechanism. 

However, we can also look at cases where the sterile 
neutrino need not be all of the dark matter, i.e., Og < 
0dm • The inferred mixing angle from the observed x-ray 
line flux is then higher, and is given by 

sin^ 26» « 7 X 10"“ (Odm/O^) . (21) 

An interesting limit to contemplate is the one where 
sterile neutrinos are produced only by the Dodelson- 
Widrow mechanism. One can then estimate the frac¬ 
tion of the total dark matter that would be constituted 
by these steriles. Solving Eqs. (20) and (21), with 
Og = Odw and 0 dm = 0.26 gives Og = Odw ~ 0.038 ~ 
0.15 0dm- Therefore, the posited 7.1 keV sterile neu¬ 
trino, even in the absence of an appreciable lepton num¬ 
ber or other nonstandard production scenarios such as 
DESNDM, could still account for roughly 15% of the 
total dark matter in this purely quantum mechanical, 
Dodelson-Widrow limit. 


D. Kinematic constraints from small- and 
large-scale structure 

As discussed in Sec. HIB, the collisionless damping 
mass/length sets the scale below which fluctuations can 
get damped by dark matter particle free-streaming. Ob¬ 
servations of large-scale structure, e.g., the Ly-a forest 
and galaxy clustering, put upper bounds on the col¬ 
lisionless damping scale. These correspond to model- 
dependent lower limits on sterile neutrino rest mass [42, 
54, 99-104]. Warm dark matter is also known to flat¬ 
ten the cores of dark matter haloes in dwarf spheroidals, 



11 


as well as decrease the expected number of low-mass 
satellites in larger dark matter haloes [105], although 
some recent work has argued that the deviation of warm 
dark matter density profiles, relative to the Navarro- 
Frenk-White (NFW) profile observed in CDM simula¬ 
tions, could be minimal [106]. Observations of stellar ve¬ 
locity dispersion profiles of dwarf galaxies can therefore 
be used to put constraints on the dark matter phase- 
space density, which again leads to mo del-dependent 
bounds on the rest mass of the candidate sterile neu¬ 
trino [107]. 

It has been suggested that observations of structure 
on small (i.e., dwarf galaxy) scales are inconsistent with 
the results of ACDM [105, 108] simulations. While it is 
not clear at this point whether the discrepancies could 
be resolved by incorporating baryonic feedback effects, 
solutions are also being sought via alterations to the 
standard CDM paradigm. For example, it has been ar¬ 
gued that a resonantly enhanced sterile neutrino, that is 
slightly warm, but nevertheless not ruled out by present- 
day observations, could lie in the sweet spot for allevi¬ 
ating some of these inconsistencies [29]. Eventually, 21- 
cm observations also may weigh in on these issues, ei¬ 
ther through flagging intergalactic medium heating from 
WIMP annihilation [109], or other insights on large-scale 
structure [110-112]. The latter indicate that perhaps 21- 
cm observations could probe structure down to length 
scales (comoving to the current epoch) of abont ~ 0.01 
Mpc, about an order of magnitude better than the scales 
probed by Lyman-a forest observations. 

Our model would likely be confronted with some of 
these constraints towards the lighter end of sterile neu¬ 
trino rest-mass range, i.e., for nis ^ a few keV, where 
the collisionless damping scales would be relevant for the 
issues discussed above. 


E. Laboratory constraints 

It is difficult to engineer direct laboratory probes of 
the candidate sterile neutrino dark matter particles and 
the candidate diluton particles we discuss here. This is 
because the couplings or rest masses of the particles in¬ 
volved can be out of reach for energies and sensitivities 
of existing or future experiments. 

Nevertheless, some laboratory probes can nibble 
around the edges of interesting parameter space for dark 
matter sterile neutrino candidates. For example, the KA- 
TRIN experiment and other direct beta decay endpoint 
experiments can target the contribution of heavy neu¬ 
trino mass eigenvalues in the coherent sum entering into 
the projection of electron flavor neutrinos in this process 
[II3-116]. 

Collider experiments, in principle, have much to say 
about beyond standard model particles and potentially 
about dark matter [117]. For example, if the dilutons 
were on the lighter end of the rest-mass range considered 
here, i.e., uih ~ TeV, then existing and future colliders 


could allow us to constrain their lifetimes [44]. Dilutons 
produced in colliders could be detected if they were to 
subsequently decay inside a detector. The detection rate 
would be proportional to the production rate times the 
ratio of crossing time to the Lorentz dilated lifetime of 
the diluton. 

Dilutons that are heavy sterile neutrinos also could be 
indirectly inferred via their impact on electroweak preci¬ 
sion observables such as the invisible Z-decay width, the 
IF-boson mass, and the charged-to-neutral current ratio 
for neutrino scattering [118]. Near-future collider experi¬ 
ments could potentially probe the effects of heavy sterile 
neutrinos on lepton-flavor-violating Z decays [119], un¬ 
fortunately their predicted sensitivities would not be high 
enough for sterile neutrinos that are sufficiently long- 
lived to be dilutons. TeV-scale sterile neutrinos could 
also influence the neutrinoless double beta decay rate 
through their contribution to the effective Majorana neu¬ 
trino mass. However, this would require a significant 
amount of fine-tuning [120], resulting in large active- 
sterile mixing and possibly rendering the sterile neutrino 
an unsuitable diluton candidate. 

Of course, current and near-future colliders are not 
likely to have the energy reach required to probe the 
physics of dilutons heavier than a few TeV in rest mass. 
High-energy cosmic ray detectors and neutrino telescopes 
such as IceCube, on the other hand, could potentially be 
useful in probing this high-energy scale physics. 


F. Compact object constraints 

A signihcant fraction of the range of rest-mass and vac¬ 
uum mixing parameters for viable sterile neutrino dark 
matter created through scattering-induced decoherence 
or resonant channels can also affect core collapse super¬ 
nova physics [3]. Conversion of active-to-sterile neutri¬ 
nos, and perhaps sometimes back again to active states, 
can affect energy and lepton number transport in the 
core, energy deposition in the mantle below the shock 
[3, 121-129], and even proto-neutron star “kicks” asso¬ 
ciated with the neutrino burst [21, 130-133]. However, 
with a small enough mixing angle, and with sufficiently 
high rest mass (i.e., well above any resonant condition in 
the core), sterile neutrino dark matter candidates consid¬ 
ered here in the DESNDM model can manage to avoid 
changing compact object physics. 


G. Differentiating between resonant production 
and DESNDM scenarios 

A relic density of sterile neutrinos comprising the 
dark matter could be produced by resonantly enhanced 
scattering-induced decoherence, as stated earlier. How¬ 
ever, for sterile neutrinos with rest masses in our range 
of interest, this would require a lepton number that is 
several orders of magnitude bigger than the baryon num- 
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ber. If future lepton number constraints, e.g., from pre¬ 
cision measurements of primordial helium and deuterium 
abundances, were to push the upper limit on the observa- 
tionally inferred lepton number to below what would be 
required for resonant production, then it would force us 
to consider alternative models, if sterile neutrinos indeed 
were to be the dark matter. Additionally, sterile neutri¬ 
nos produced resonantly would have warmer energy spec¬ 
tra, i.e., larger collisionless damping scales, compared to 
dilution generated sterile neutrinos of the same mass. 
Improved constraints on the collisionless damping scale 
from future Ly-a observations and possibly 21-cm obser¬ 
vations, as well as an improved understanding of small- 
and large-scale structure formation through sophisticated 
simulations, could lead to certain models gaining favor 
over others. 


V. CONCLUSIONS 

We have described a generic mechanism, the DESNDM 
model, whereby sterile neutrinos with rest masses in the 
range ~ keV to ~ MeV could acquire relic densities that 
allow them to be the dark matter. The key assumptions 
of this model are that: (1) the dark matter candidate 
sterile neutrino is in thermal and chemical equilibrium at 
very high temperature scale in the early Universe, (2) this 
particle decouples at very early epochs, (3) subsequent 
to this dark matter candidate sterile neutrino decoupling 
event there is prodigious entropy generation from the out- 
of-equilibrium decay of a different particle, the diluton, 
and (4) the diluton is presumed to be very massive and 
to possess nonrelativistic kinematics when it decays, but 
is also assumed to have been in thermal and chemical 
equilibrium at very early epochs, possessing relativistic 
kinematics at the time of its decoupling. 

We do not identify the diluton with a specific parti¬ 
cle candidate. Instead, we consider the generic issues in¬ 
volved with out-of-equilibrium particle decay in the early 
Universe and attendant cosmological, observational, and 
laboratory constraints. This leads us to consider dilu- 
tons with rest masses in the ~ TeV to ~ EeV range, 
possessing rather long decay lifetimes. Candidate dilu¬ 
ton particles might include, for example, heavy sterile 
neutrinos, different from the ones we might consider for 
the dark matter particle, and supersymmetric particles 
with R-parity-violating decays into standard model par¬ 
ticles. The latter scenario would suggest a very high su¬ 
persymmetry scale, and a novel and heterodox role for 
supersymmetric particles in the dark matter problem. 

The DESNDM dilution mechanism for producing the 
sterile neutrino relic densities of interest necessarily re¬ 
sults in a corresponding relic energy spectrum which, 
though thermal in shape, can be quite cold compared 
to a standard energy spectrum characterized by a relic 
photon or relic active neutrino temperature. This allows 
our sterile neutrino dark matter to behave like CDM in 
many cases, even though the actual rest masses of the 


sterile neutrinos are modest. Additionally, the sterile 
neutrino population could acquire a “warm,” albeit sub¬ 
dominant component arising via scattering-induced de¬ 
coherence, even in the zero lepton number limit, which 
could have some observable effects on structure formation 
on the small scales. 

Interestingly, the DESNDM model for generating ster¬ 
ile neutrino relic densities can be, effectively, nearly inde¬ 
pendent of the vacuum mixing angle characterizing the 
mixing of this sterile neutrino with any of the active neu¬ 
trino flavors. The model requires only that this mixing 
angle be smaller than that required to effect population 
of a sterile neutrino sea from the seas of active neutrinos 
in the very early Universe. This is unlike other models 
for producing a sterile neutrino relic dark matter den¬ 
sity. For example, in scattering-induced decoherence and 
resonant enhancement of this process, the vacuum mix¬ 
ing angle is a key parameter, so that sterile neutrino rest 
mass, this mixing angle, and perhaps other parameters 
like lepton number, uniquely determine the relic density. 
This means that x-ray observational constraints can, in 
principle, definitively rule out ranges of sterile neutrino 
rest mass and vacuum mixing parameter space. 

In the DESNDM model the relic density is set by differ¬ 
ent physics. The existence of a nonzero vacuum mixing 
with active neutrinos will, of course, still guarantee a ra¬ 
diative decay channel for this sterile neutrino particle. 
However, in the DESNDM model a given sterile neutrino 
rest mass with a dark matter relic density need not have 
a mixing angle large enough to produce an x-ray flux 
sufficient for detection. As a consequence, the DESNDM 
mechanism can evade all x-ray bounds, and the detection 
of a dark matter sterile neutrino decay line would be a 
lucky, but not inevitable development. 

That does not mean that there are no potential obser¬ 
vational or experimental handles on sterile neutrino dark 
matter produced via the DESNDM mechanism. First, 
note that dilution can make dark matter sterile neutri¬ 
nos with a wider range of rest masses than is possible 
in scattering-induced decoherence models. As discussed 
above, new experiments like NuSTAR can probe higher- 
energy X rays. At x-ray energies above ^ 10 keV, the 
expected x-ray backgrounds are lower than they are in 
the “sweet spot” of a few to 10 keV for the XMM and 
Chandra x-ray telescopes. Second, as discussed above, 
large-scale structure simulations and observations may 
be able to produce finer probes of the dark matter char¬ 
acter. 
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